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ABSTRACT 

We investigate the power of geometrical estimators on detecting non-Gaussianity in 
the cosmic microwave background. In particular the number, eccentricity and Gaussian 
curvature of excursion sets above (and below) a threshold are studied. We compare 
their different performance when applied to non-Gaussian simulated maps of small 
patches of the sky, which take into account the angular resolution and instrumental 
noise of the Planck satellite. These non-Gaussian simulations are obtained as per- 
turbations of a Gaussian field in two different ways which introduce a small level of 
skewness or kurtosis in the distribution. A comparison with a classical estimator, the 
genus, is also shown. We find that the Gaussian curvature is the best of our estima- 
tors in all the considered cases. Therefore we propose the use of this quantity as a 
particularly useful test to look for non-Gaussianity in the CMB. 

Key words: methods: statistical - cosmic microwave background 



1 INTRODUCTION 



Future cosmic microwave background (CMB) experiments, 
such as the VSA (Jones & Scott 1998), the NASA MAP 
. !^ ' satellite (Bennett et al. 1997) and the Planck mission from 
K^ , ESA (Tauber 2000), will provide with unique data to con- 
5_j straint fundamental cosmological parameters as well as to 
Cu distinguish between competing theories of structure forma- 
tion in the early universe. In particular, the standard infla- 
tionary model predicts Gaussian fluctuations whereas topo- 
logical defects give rise to non-Gaussian signatures in the 
CMB. Therefore, the study of the Gaussianity of the CMB 
is a key issue to understand the nature of the primordial per- 
turbations that led to the formation of large-scale structure. 
In addition, future high resolution maps should be searched 
for any traces of non-Gaussianity since both foregrounds and 
systematic errors can also leave non-Gaussian imprints. 

A large number of methods have been proposed in the 
literature to test the Gaussianity of the CMB (for a review 
see Barreiro 2000). They include extrema correlation func- 
tion (Kogut et al. 1996, Barreiro et al. 1998, Heavens & 
Sheth 1999), properties of hot and cold spots (Coles & Bar- 
row 1987, Martinez-Gonzalez et al. 2000), Minkowski func- 
tionals (Coles 1988, Gott et al. 1990, Kogut et al 1996), 
bispectrum analysis (Ferreira et al. 1998, Heavens 1998, 
Magueijo 2000), wavelets (Pando et al. 1998, Hobson et al 
1999, Mukherjee et al. 2000, Barreiro et al. 2000), muhifrac- 



tals (Pompilio et al. 1995) and partition function (Diego et 
al. 1999, Martinez-Gonzalez et al. 2000). 

In this paper we investigate the use of geometrical es- 
timators on detecting non-Gaussianity in the CMB using 
small scale simulated maps (for Gaussian and non-Gaussian 
flelds). These simulations take into account the angular res- 
olution and sensitivity predicted for two of the Planck fre- 
quency channels (the LFI 100 GHz and the HFI 217 GHz 
channels). In particular, we study the performance of the 
number, eccentricity and Gaussian curvature of excursion 
sets above (or below) a threshold. These quantities have 
been studied for a Gaussian field by Barreiro et al. (1997). 
For comparison, we also present the same analysis for the 
genus, a topological quantity that has been proposed as a 
good estimator to test the Gaussianity of the CMB. 

Simulated maps for the standard inflationary model can 
be very accurately (at the level of ~ 1 per cent in the power 
spectrum) obtained (e.g. Seljak & Zaldarriaga 1996, Hu et 
al. 1998). However this is not the case for models based on 
topological defects. Besides, the power spectrum predicted 
by pure topological defects scenarios (see e.g. Pen et al. 1997, 
Contaldi et al. 1999) seems not to be in agreement with the 
latest CMB data produced by BOOMERANG (de Bernardis 
et al. 2000) and MAXIMA (Hanany et al. 2000) but the pos- 
sibility of an hybrid scenario that combines both inflation 
and topological defects remains viable (Bouchet et al. 2000, 
Contaldi 2000) . To take into account for all of these uncer- 
tainties, we have used two generic non-Gaussian flelds. In 



© 0000 RAS 



R.B.Barreiro et al. 



particular, physically motivated non-Gaussian models, e.g. 
those based on topological defects or non-standard infla- 
tion, generically produce certain level of skewness or/and 
kurtosis. Since these are the lowest order kumulants they 
are also the least affected by the presence of instrumental 
white noise. Therefore, we have generated two different kind 
of non-Gaussian maps obtained as weak perturbations of a 
Gaussian field. The first one, that we call P distribution, 
has an asymmetric 1-point density function (1-pdf) with a 
positive tail, what introduces a small level of skewness. The 
1-pdf of the second non-Gaussian field, the B distribution, is 
symmetric and have a higher peak and longer tails than the 
Gaussian one. This produces a positive value of the (excess) 
kurtosis. Some of these characteristics have been shown to be 
produced by topological defect models by different authors. 
Turok (1996) performed sub-degree scale simulations of the 
CMB from cosmic defects, finding a clear excess of hot spots 
for monopoles and textures, a signature also present on our 
P distribution. Perivolaropoulos (1993) used an analytical 
model to investigate the statistics of the CMB temperature 
maps induced by topological defects. He obtained a 1-pdf for 
cosmic strings that presents a central peak and long tails, 
resembling our B distribution. 

The outline of this paper is as follows. In §2 we discuss 
the geometrical quantities that we test as non-Gaussian es- 
timators. §3 describes our Gaussian and non-Gaussian CMB 
simulations. The results, including the effect of instrumental 
noise, are presented in §4. Finally, the main conclusions of 
our work are summarised in 55. 



2 GEOMETRICAL ESTIMATORS 

We have already pointed out the importance of testing the 
Gaussianity of the CMB. An interesting possibility is the use 
of geometrical estimators related to excursion sets above or 
below a given threshold. In particular, we have considered 
the mean number, eccentricity and Gaussian curvature of 
regions above or below a threshold u (in units of the signal 
dispersion). All these quantities can be easily estimated from 
a map. The procedure is as follows. First the excursion sets, 
i.e. connected pixels, above v (and below —u) are identified. 
Each of them is fitted to a paraboloid with parameters a, b 
and c (see Fig. hi). The eccentricity e and Gaussian curvature 
K, of each of the spots are then estimated as 



1- 



ba 



(1) 



where a, b are in arcminutes and c in uints of the map dis- 
persion. Only those excursion sets with a single maximum 
(or minimum when below the threshold) and formed by at 
least eight pixels are considered to obtain the mean Gaus- 
sian curvature and eccentricity of the map. 

In addition, these quantities can be ob- 
tained (semi)analytically for a homogeneous and isotropic 
Gaussian field, at least locally around the maximum. On 
the one hand, the probability density functions of the num- 
ber, eccentricity and Gaussian curvature of maxima (which 
coincide asymptotically with the excursion sets) for a ho- 
mogeneous and isotropic Gaussian field have been studied 




Figure 1. To obtain the mean eccentricity and Gaussian cur- 
vature of the excursion sets of the CMB, these are fitted to a 
paraboloid as the one shown in the figure, a and b are the major 
and minor axes of the ellipse that forms the basis, c is the distance 
between the geometrical centre of the ellipse and the maximum 
of the paraboloid 



by Barreiro et al. (1997). On the other hand, the expected 
number of excursion sets < N > over the whole sphere can 
be estimated as (Vanmarcke 1983): 



<iV> = 



■^0c eric{u/V2) 



(2) 



where 6c is the coherence angle of the field that depends 
only on the 2-point correlation function. 

For comparison, we have also extended our analysis to 
the genus, a topological quantity that has been commonly 
used to test the Gaussianity of the CMB. The genus at a 
given threshold can be approximately estimated as the num- 
ber of isolated high-temperature regions minus the number 
of isolated low-temperature regions. The mean genus < g > 
per unit area for a homogeneous and isotropic Gaussian field 
is given by 



<g>= 



(27r)3/2ei 



-1.^/2 



(3) 



To estimate this quantity from a map we have used a method 
based on adding up the contribution to the genus from each 
vertex of the pixelised map. This method is described in 
detail in Jones (1999) (see also Hamilton et al., 1986). 

In order to study the validity of our algorithms on es- 
timating the above quantities, we have compared the val- 
ues computed from simulations with the theoretical expec- 
tations for the Gaussian case. Fig. H shows the mean number 
of regions, eccentricty, the inverse of the Gaussian curvature 
and genus obtained from 500 Gaussian realisations for a A- 
CDM model. The size of the maps is 512*512 pixels, with a 
pixel size of l'.5 and they have been smoothed with a Gaus- 
sian beam of FWHM=5.'5 . The solid lines correspond to 
an approximation of the theoretical expectations for each 
of these quantities. For the number of regions and genus, 
the curves have been obtained from equations ti and H re- 
spectively, whereas the inverse of the Gaussian curvature 
and the eccentricity were calculated using the probability 
density functions given by Barreiro et al.(1997). We can see 
that there is a certain difference between the computed and 
predicted values for the number of regions, especially at low 
thresholds. This is due to the fact that the Vanmarcke con- 
jecture is an approximation that overestimates the number 
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Figure 2. Comparison of the average values of our estimators 
as computed from 500 simulated maps (open circles) with an ap- 
proximation of the theoretical expectations for a Gaussian field 
(solid lines). The error bars show the dispersion obtained from 
the simulations. The panels correspond to (a) number of regions, 
(b) eccentricy and (c) inverse of the Gaussian curvature for ex- 
cursion sets above (or below for negative v) a threshold v, and 
(d) genus. 



of excursion sets at low thresholds. Regarding the inverse of 
the Gaussian curvature as well as the eccentricity, we have 
to point out that the theoretical values have been obtained 
locally around the maxima, whereas on the simulations we 
consider the whole of the excursion set and discard those 
with more than a single maximum and less than eight pixels. 
This procedure biases the expected values and can explain 
the discrepancies between the theoretical and computed val- 
ues, that are however reasonably close. Finally, as expected, 
the values computed from the simulations for the genus are 
in good agreement with the theoretical curve. 



3 SIMULATIONS 

In order to test the performance of our estimators, we have 
simulated small patches of the sky of size 12.8 square de- 
grees with Gaussian and non-Gaussian statistical properties. 
We have considered two different angular resolutions that 
correspond approximately to two of the Planck frequency 
channels. On the one hand, we have produced maps with 
pixel size of 3 arcminutes, smoothed with a Gaussian beam 
of FWHM=10', which corresponds to the resolution of the 
LFI 100 GHz channel. On the other hand, we have simulated 
maps with the characteristics of the HFI 217 GHz chan- 
nel, i.e., pixel size of 1.5 arcminutes and Gaussian beam of 
FWHM=5'.5. We have also taken into account the instru- 
mental noise of these channels. 

The Gaussian simulated maps have been generated from 



a CDM model normalised to COBE/DMR and parameters 
VLm = 0.3, Q.K = 0.7, h = 0.5, fib = 0.05 and n = 1. The 
power spectrum of such a model has been obtained with 
CMBfast (Seljak & Zaldarriaga 1996). The non-Gaussian 
simulations have been produced as follows. First, an uncor- 
related 2D Gaussian map with zero mean and unit variance 
is obtained. This map is then transformed in two different 
ways according to the following expressions (Weinberg & 
Cole 1992): 



Fi 

F2 



^(e"^ -t- e-°^) 
2[l-|-(l + 4a2)e2°^] 



P transformation (4) 
B transformation (5) 



where G corresponds to the initial Gaussian field and F\ ,F2 
to the transformed fields. The level of non-Gaussianity in- 
troduced can be controlled with the parameter a. In partic- 
ular when a ^ we recover the Gaussian case. In this way, 
we obtain two kind of uncorrelated non-Gaussian fields with 
zero mean and unit variance. Next, we rescale the amplitude 
of each Fourier mode to obtain the same power spectrum as 
for the Gaussian CDM model. By rescaling the power spec- 
trum the 1-pdf is modified with respect to the uncorrelated 
one, but it is still non-Gaussian (hereinafter, when we re- 
fer to P and B transformations, it should be understood 
the transformations given by (4) and (5) plus a subsequent 
rescaling of the power spectrum as explained before). In par- 
ticular, for the P transformation we are constructing a non- 
Gaussian field whose distribution is skewed with a positive 
tail, whereas for the B transformation, the distribution is 
broadened with both positive and negative tails. Finally, the 
new maps are smoothed with the required Gaussian beam. 
For both kind of non-Gaussian fields, we have considered 
small enough values of a such that their 1-point distribu- 
tions are indistinguishable from the Gaussian one. In addi- 
tion, the Gaussian and non-Gaussian simulations share the 
same power spectrum by construction. Fig. |3| (top panels) 
shows the 1-pdf for the Gaussian and non-Gaussian fields for 
a = 0.8 that have been obtained averaging over 500 simula- 
tions with angular resolution corresponding to the 217 GHz 
channel. Noise has not been included in this case. For the 
sake of clarity, we have not plotted the value of the disper- 
sion at each bin of the histogram as error bars but they com- 
pletely overlap for the Gaussian and non-Gaussian 1-pdf's. 
We can see that the P distribution has a positive skewness, 
whereas the B one presents a small level of kurtosis. This 
can also be seen in the middle and bottom panels, where the 
histogram of the values of the skewness and kurtosis have 
been plotted for the 500 Gaussian and non-Gaussian simu- 
lations (these quantities have been obtained using unbiased 
estimators as explained in Hobson et al. 1999). Table hi gives 
the average values as well as the dispersions of the skewness 
and kurtosis distributions obtained from the previous sim- 
ulations. From these values, we can see that the Gaussian 
and non- Gaussian distributions can be differentiated only at 
the ~ Icr level using the skewness and kurtosis. 
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Figure 3. Top: 1-pdf function for tlie Gaussian (solid line), P 
(dashed line) and B (dotted line) distributions (with unit dis- 
persion) obtained averaging over 500 simulations. The parame- 
ter a = 0.8 has been used for the transformed fields. The size 
of the simulated maps was 12.8 square degrees with a pixel of 
1.5 arcminutes and they were smoothed with a Gaussian beam 
with FWHM=5.5 arcminutes, which corresponds to the angular 
resolution of the Planck 217 GHz channel. The corresponding his- 
tograms for the skewness (middle) and kurtosis (bottom) for the 
same simulations are also plotted. 



Table 1. Values of the mean and dispersion of the skewness and 
kurtosis distributions obtained from 500 simulations. 



Distribution 



<53> 



<^(S3) 



<Ki> 



criKi) 



G 

P (a=0.8) 
B (a=0.8) 



3.72 X 10" 

0.188 
7.35 xlO" 



7.80x10" 
8.60 xlO 
9.30 xlO 



-6.55 XlO' 
0.117 
0.270 



0.125 
0.161 
0.358 



4.1 Noiseless case 

First we will consider simulated maps where no noise has 
been included. The top panels of Fig. Ek show the mean 
value and dispersion of the number of spots above (for 
positive v) or below (for negative i/) a threshold obtained 
over 500 Gaussian and non-Gaussian realisations of size 
256 X 256 (pixel 3') and smoothed with a Gaussian beam 
of FWHM=10' (LFI 100 GHz channel). Open circles, sohd 
triangles and solid squares correspond to the Gaussian, P 
and B fields respectively. In order to establish a better com- 
parison we have constructed the likelihood function for the 
Gaussian and non-Gaussian maps (compared to the Gaus- 
sian model) given by 



-i(a;=-<a:"3>)^M-l(a; 



2tt\M\ 



1/2 



M,, = 



<{xf-<xf >)ixf- 



< X]" >) > 



where s refers to the Gaussian and non-Gaussian simulated 
data and G to the Gaussian ones. M is computed from the 
Gaussian realisations. The histogram of the likelihood val- 
ues for our Gaussian and non-Gaussian fields when using 
the number of spots are plotted in the bottom panels of 
Fig. Ha. The mean values, dispersions and likelihood curves 
for the eccentricity, Gaussian curvature and genus (in abso- 
lute value) are given in Figs. Wb, Mr andWd respectively. For 
negative thresholds the Gaussian curvature and eccentricity 
correspond to the excursion sets below that threshold. We 
can see in these figures that the P distribution can be clearly 
distinguished from the Gaussian model using the number, 
Gaussian curvature or the genus since the likelihood curves 
do not overlap, whereas the eccentricity does not discrimi- 
nate between both distributions. In the three first cases, the 
100% of the non-Gaussian simulations lie outside the 99% 
confidence limit of the Gaussian case, whereas only 23.8% of 
the P simulations could be distinguished from the Gaussian 
model at the 99% c.l. using the eccentricity. To establish a 
better comparison in the performance of our estimators, we 
' have calculated a factor d defined as 



d^ 



ma — mjvG 



(7) 



4 RESULTS 

In order to compare the power of our estimators, we have 
obtained the mean number, eccentricity, Gaussian curvature 
and genus for a large number of Gaussian and non-Gaussian 
simulations with two different angular resolutions which cor- 
respond to the LFI 100 GHz and the HFI 217 GHz channels 
of the Planck satellite. We have chosen a value of a = 0.8 
for the non-Gaussian maps, that introduces small but in- 
teresting differences. We have also investigated the effect 
of instrumental noise. In order to have a reasonable num- 
ber of spots at any given threshold in all of our realisations 
we have considered thresholds in the range —2.5 < v < 2.5 
(—2 < u < 2.5 for the P distribution in some cases). In most 
of our cases thresholds close to i' = do not contain very 
relevant information and have not been included to speed 
up the numerical calculations. 



where mc and mjvG correspond to the median of the likeli- 
. hood curve of the Gaussian and non-Gaussian cases respec- 
tively and ac is the dispersion of the likelihood curve for the 
Gaussian simulations. Therefore this quantity is an indica- 
tion of how well we can separate between the Gaussian and 
non-Gaussian distribution with a given estimator. Table ti 
gives the values of this quantity for the different cases con- 
sidered in this paper without including instrumental noise. 
In the same table, the quantity / indicates the percentage 
of non- Gaussian simulations that lies outside the Gaussian 
99% c.l. for a given estimator. We can see that for the P 
distribution in the case of the LFI 100 GHz channel, the 
highest value of d (and therefore the best estimator) corre- 
sponds to the Gaussian curvature (d=44.8 versus d=24.7 
for the genus, the next best quantity). Regarding the B 
model, the likelihood curves overlap for all the estimators, 
except for the Gaussian curvature, that gives again the best 
results (d=26.1, /=100%). The other three estimators can 
not discriminate between both distributions, being the ec- 
centricity (d=3.9, /=60.6%) slightly better than the genus 
(d=3.4,/=56.2%) and the number (d=2.6,/=45.8%). 
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Figure 4. A comparison of the power of our estimators on detecting small scale non-Gaussianity using simulated maps with the 
specifications of the Planck 100 GHz channel (no noise) is shown. The four figures in the top-left corner (labelled (a)) correspond to 
the results for the number of spots above (for positive thresholds) or below (for negative thresholds) u. The two top figures of this 
set show the mean values and dispersions obtained from the distribution of the number of excursion sets over 500 simulations for the 
Gaussian (open circles), P (solid triangles) and B (solid squares) fields. The histogram of likelihood values computed from the same 
simulations are plotted in the two corresponding bottom figures for the Gaussian (solid line), P (long-dashed line) and B (short-dashed 
line) distributions. The same results are given for the eccentricity (top-right figures, labelled (b)), Gaussian curvature (lower-left figures, 
(c)) and genus (lower-right figures, (d)). 



The results for the small simulated maps of the 217 
GHz channel (512 x 512 maps, pixel=l'.5, FWHM=5'.5) 
are shown in Figs, ga (number), ^b (eccentricity), |5|c (cur- 
vature) and Hd (genus). As before, the P distribution can 



be clearly distinguished from the Gaussian one using the 
Gaussian curvature, genus or number of excursion sets but 
not with the eccentricity. In particular the Gaussian cur- 
vature is again the best of our statistics with d=67.1 and 
f—100% whereas the mean number (d=28.5, /=100%) and 
genus (d=24.1, f—100%) produce very similar results. The 
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Table 2. Values of d and / (no noise) 
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Figure 5. As in Fig.l4 but with the results corresponding to the 217 GHz channel. 
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< K > also succeeds in distinguishing between the B and 
Gaussian distributions (d=24.9,/=100%), faihng the rest of 
estimators to discriminate between them. 

It is interesting to point out that, in spite of the better 
resolution of the 217 GHz channel, the level of success of 
our estimators is very similar for both the 100 and 217 GHz 
channels. This can be understood taking into account that 
the coherence angle for the considered A- model is ~ 10'. 
Therefore the information about the intrinsic structure of 
the CMB contained in both channels is very similar. 

4.2 The effect of instrumental noise 

We have also investigated the effect of instrumental noise 
in the performance of our estimators. In particular, we have 
included white Gaussian noise at the level predicted for the 
Planck satellite after 12 months of observation for the two 
channels that we are considering (AT/T per resolution el- 
ement = 4.3 X 10~^ in both cases). For our model this 
corresponds to a signal to noise ratio (S/N) per pixel of 
S/N ~ 2.8 and S/N ~ 2.7 for the 100 and 217 GHz chan- 
nels respectively. In order to increase the S/N, the simulated 
signal plus noise maps have been smoothed with a Gaussian 
beam with FWHM equal to the antenna of the correspond- 
ing channel. 



The results are summarised in Table and Figs. Ma (Gaus- 
sian curvature, 100 GHz), 0b (genus, 100 GHz), He (Gaus- 
sian curvature, 217 GHz) and nd (genus, 217 GHz). Only the 
two best estimators are shown. The most interesting point 
is that again the Gaussian curvature is the best estimator 
in all the considered cases. In particular for the P distribu- 
tion the 100% of the realisations have lied outside the 99% 
confidence level of the Gaussian model for both frequency 
channels. The genus and the number of excursion sets give 
again very similar results. They are able to discriminate be- 
tween most of the P and Gaussian realisations, although the 
confusion level is larger than for the Gaussian curvature. The 
eccentricity completely fails to distinguish between them. 

For the B distribution, there is a small overlap be- 
tween both likelihood curves when using the Gaussian cur- 
vature but still most of the realisations fall outside the 
Gaussian 99% c.l. (d=6.8,/=93% for the 100 GHz chan- 
nel, d=8.2,/=89% for the 217 GHz channel). The number, 
genus and eccentricity can not discriminate between the B 
and Gaussian distributions having values of d ~ 1 . 



5 CONCLUSIONS 

We have performed a comparative study of the ability of 
three geometrical estimators (mean number, eccentricity 
and Gaussian curvature of excursion sets above or below 
a threshold) as well as the genus on detecting small scale 
non-Gaussianity in the CMB. To test such estimators, we 
have generated two kind of non-Gaussian realisations that 
are perturbations of a Gaussian field. The first one (P dis- 
tribution) introduces a certain level of skewness in the 1-pdf 
whereas the second one {B distribution) has a non-negligible 
kurtosis. This could resemble the characteristics of some 
topological defect models. The simulated maps are small 



patches of the sky (12.8 square degrees) with angular reso- 
lutions which correspond approximately to two of the Planck 
satellite frequency channels, the LFI 100 GHz channel (pixel 
size of 3', FWHM=10') and the HFI 217 GHz (pixel size of 
1'.5, FWHM=5'.5). We have also taken into account the 
expected level of instrumental noise in these channels. 

We find that the Gaussian curvature is the best of our 
estimators. This is a robust conclusion, since this quantity 
performs considerably better than the genus, number and 
eccentricity of excursion sets for both types of non-Gaussian 
fields, in presence or absence of noise and for the two consid- 
ered angular resolutions. Regarding the number of excursion 
sets and the genus, these two quantities carry very similar 
information, as one would expect. They can discriminate be- 
tween the P and Gaussian models at a fairly high level (al- 
though always worse than the Gaussian curvature) in all the 
considered cases. However, these quantities fail to discrim- 
inate between the B and Gaussian distributions, being the 
level of confusion especially high when instrumental noise is 
introduced. The performance of the eccentricity is very poor. 
In none of the considered cases can clearly discriminate be- 
tween the non-Gaussian and Gaussian fields, although it is 
interesting to point out that it does perform at similar or 
better levels than the genus and the number of excursion 
sets when testing the B distribution. Therefore, the main 
conclusion of our work is that the mean Gaussian curvature 
of excursion sets could be a very useful estimator to look for 
small scale non-Gaussianity in the CMB. 
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Figure 6. The results for the Gaussian curvature-100 GHz (four top-left figures, (a)), genus-100 GHz (top-right figures, (b)), Gaussian 
curvature-217 GHz (bottom-left figures, (c)) and genus-217 GHz (bottom-right, (d)) are shown. Noise has been included in the simulations. 
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